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An observer surrounded by sufficiently small spherical light sources at a fixed dis- 
tance will see a pattern of elliptical images distributed over the sky, owing to the 
distortion effect (shearing effect) of the spacetime geometry upon light bundles. In 
lowest non-trivial order with respect to the distance, this pattern is completely deter- 
mined by the conformal curvature tensor (Weyl tensor) at the observation event. In 
. this paper we derive formulae that allow one to calculate these distortion patterns in 

terms of the Newman-Penrose formalism. Then we represent the distortion patterns 
graphically for all Petrov types, and we discuss their dependence on the velocity of 
the observer. 

o " 
o 

O , Suggested short title: Classification of image distortions 

cr 

99; 1 Introduction 

> ■ 

^ ■ The general-relativistic light deflection, i.e., the influence of the spacetime geometry on 
the paths of light rays, has the effect that, in general, the apparent shape of a distant 
object at the celestial sphere of an observer will be distorted. This distortion effect is 
directly related to the shear of light bundles, whereas the apparent size and the apparent 
brightness of images are related to the expansion of light bundles. By the well-known Sachs 
equations Q, the shear is an effect of the conformal curvature tensor (Weyl tensor), whereas 
the expansion is influenced by the Ricci tensor along the light rays. In the following we 
concentrate on the distortion effect, i.e. on the shear. In a pioneering paper, Kristian and 
Sachs [fj] introduced a certain measure for the distortion effect, given in terms of a series 
expansion with respect to the distance between source and observer, and suggested using 
it as a cosmological observable. In lowest non-trivial order with respect to the distance the 
distortion is completely determined by the Weyl tensor at the observation event. It is then 
a natural question to ask to what extent the distortion is already determined by the Petrov 
type of the Weyl tensor. Up to now only some special aspects of this question have been 
addressed in the literature (see Penrose and Rindler ||, Volume II, Chapter 8); it is the 
purpose of this paper to give a more detailed discussion, including graphical representations 
of the distortion patterns on the sky for each Petrov type. 

To gain understanding the distortion effect it is didactically helpful to begin the discus- 
sion with an over-idealized universe in which all galaxies are perfectly spherical. We shall 
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assume that all light rays that reach the observer from a particular galaxy can be viewed as 
'infinitesimally close to each other', i.e., that the laws of light propagation can be linearized 
around a central light ray. Then the image of each galaxy results by applying a linear map 
to a circle, i.e., it is an ellipse. By the term 'distortion pattern' we mean the distribution 
of those ellipses over the celestial sphere of the observer, for galaxies at a specific distance. 
(We are interested only in the shape, not in the size, of the ellipses.) Having fixed the 
distance, the distortion pattern is given by a non-negative scalar function on the celestial 
sphere that gives the eccentricity of the respective ellipse, and by a direction field on the 
celestial sphere that indicates the major axis of the respective ellipse. The relevant formulae 
for determining these fields, which are due to Kristian and Sachs , will be given in Section 
H below. We shall restrict our discussion to the case that only terms of lowest non-trivial 
order with respect to the distance are to be taken into account. (In this approximation it 
does not matter which of the various non-equivalent general-relativistic notions of 'distance' 
we use, see Section below.) The distortion pattern is then completely determined by the 
Weyl tensor at the observation event, with the eccentricities depending quadratically on 
the distance and the direction field being independent of the distance. If the Weyl tensor 
is zero, then there is no distortion effect, i.e., the eccentricity is everywhere zero and the 
direction field is undetermined. If the Weyl tensor is not zero, then the distortion effect 
vanishes at exactly four (not necessarily distinct) points on the sky, corresponding to the 
four principal null directions of the Weyl tensor. Kristian and Sachs concentrate on the 
magnitude of the distortion effect, i.e., on the eccentricities which they measure in terms of 
a function they denote by e. Penrose and Rindler ||, Volume II, Chapter 8, on the other 
hand, concentrate on the direction field which they call the 'fingerprint' of the Weyl tensor. 
It is our goal to discuss both the magnitude and the direction of the distortion effect; there- 
fore, we assign to each point of the observer's celestial sphere a 'distortion length element' 
whose length gives the eccentricity and whose direction gives the major axis of the respective 
ellipse. In Section ^| we show how these quantities can be calculated with the help of the 
Newman-Penrose formalism. In Section |] we discuss, in terms of pictures, the 'distortion 
length element field' and its dependence on the oberserver's velocity for each Petrov type. 

If we want to link these deliberations with the real universe we have to face the problem 
that galaxies are not spherical. We may restrict to those galaxies which appear elliptical 
on the sky, to within reasonable approximation, and we may assume that the actual shape 
of those galaxies is approximately that of a rotational ellipsoid. The problem is that we 
do not know the actual eccentricity and the actual direction of the rotation axis for any 
individual galaxy. As long as this is true, it is in principle impossible to measure the 
distortion effect upon any individual galaxy. The only way out of this difficulty is to use 
statistical methods in order to measure the distortion effect upon sufficiently large samples 
of galaxies. This approach is based on the assumption that, on a sufficiently large scale, 
the rotation axes of galaxies are randomly distributed in the universe. (This is a reasonable 
working hypothesis, but it is not beyond any doubt. In a rotating universe, e.g., the angular 
momenta of galaxies are expected to be preferably aligned with the universal rotation axis.) 
The distortion pattern, as introduced above, should then be obtained by averaging galaxy 
ellipticities over sufficiently large fields in the sky. In this statistical sense, the distortion 
effect has actually been measured by several groups. In the following paragraphs we give a 
brief overview of the observational situation. 

The first attempt of observing the distortion effect was done already in 1966 by Kristian 
PJ with a small number of galaxies in several clusters. The observation was unsuccessful 
which, together with some reasonable assumptions, yielded an upper limit for the magnetic 
part of the Weyl tensor. This null result was confirmed in 1983 by Valdes et al. H with 
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a 300fold improvement in accuracy over Kristian, owing to a bigger number of galaxies 
and including galaxies at larger distances. The first successful observation of the distortion 
effect by statistical methods was reported in 1990 by Tyson et al. || who considered the 
distorting effect of rich clusters of galaxies upon background galaxies. The same mechanism 
- background galaxies distorted by the gravitational field of an intervening cluster - is 
usually considered as an explanation for the socalled giant luminous arcs which have been 
discovered in great number since 1986, beginning with Lynds and Petrosian [0] and Soucail 
et al. f|. (In the latter case, the linear approximation is, of course, not applicable, i.e., it is 
not justified to view all light rays coming from the source to the observer as 'infinitesimally 
close' to some central light ray.) In this sense, the distortion effect produced by galaxy 
clusters was a well-established phenomenon by the mid-1990s. On the other hand, by this 
time all attempts to observe the distortion effect produced by large-scale structure had 
failed, see Mould et al. for an unsuccessful attempt in 1994. 

With the recent advancement of CCD mosaic cameras the observational situation has 
greatly changed. Now it is possible to observe 'great' parts (~ 30 x 30 square arc-minutes) 
of the sky with a large number of galaxies (~ 10 5 ) in a single picture and to handle the 
contained information automatically. The unwrought data can be cleaned by several tech- 
niques described in HTT| , [T2| , |i~3| , [14}| which improves the significance of the effect. The results 
of these recent observations can be summarized in the following way. A non-zero distor- 
tion effect on a large scale has been measured by demonstrating a correlation of galaxy 
ellipticities in several selected fields in the sky [pi], [TT], [12|, 13|, |i"4|, and there is evidence 
that this correlation decreases if the field size is increased O, [T2L O, ITJJ . Image distortion 



by large-scale structure is often called cosmic shear, refering to the shear induced in light 
bundles on a cosmic scale. This term should not be confused with a hypothetical shear of 
the Hubble flow which occasionally is also called 'cosmic shear' or 'cosmo logical shear'. 

For future perspectives of observations we refer to the Sloane Digital Sky Survey (see 
tittp : / /www . sdss . org| ) and to the Deep Lens Survey (see [http : / /dls . bell-labs . com| ) 
which are both under way, but also to the proposed Dark Matter Telescope (see, e.g., 
Tyson, Wittman and Angel |15|]). It does not seem to be overly optimistic to hope that in 
some years the observational material on the distortion effect covers large parts of the sky. 

Observations of image distortions are usually evaluated with the help of the weak-lensing 
formalism which is based on early ideas of Gunn 
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(also see Webster [18[]) and was 

van 



developed since 1990 by several authors including, e.g., Miralda-Escude [[19], Kaiser pO 
Waerbeke 



Bernardeau and Mellier 
and Schneider Wl 
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for a comprehensive review we refer to Bartelmann 
This formalism uses Newtonian approximations in a weakly perturbed 
Friedmann-Robertson- Walker model. As the Weyl tensor of a Friedmann-Robertson- Walker 
model vanishes, the distortion effect in such a universe is produced by the perturbations 
alone. Among other things, the weak-lensing formalism provides a relation between the 
correlation function of ellipticities and the power spectrum of mass fluctuations, thereby 
relating the distortion effect to the distribution of (dark) matter. These results rely, of 
course, on the Newtonian approximations. The fact that the distortion effect is determined 
by the Weyl tensor, on the other hand, is purely geometric and quite general. In the weak- 
lensing formalism this central role of the Weyl tensor is somewhat disguised and partly 
overshadowed by the approximations. In particular, the effect of the magnetic part of the 
Weyl tensor is completely neglected because it has no Newtonian counterpart. The present 
paper is motivated, at least partly, by our desire to rediscuss the role of the Weyl tensor 
for image distortions in view of the recent observations. More precisely, we want to ask if 
the observed image distortion by large-scale structure can be used to gain some information 
about the Weyl tensor of the universe. We shall discuss this question, on the basis of the 
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results derived in the body of this paper, in the conclusions. 

Leaving aside all applications to cosmology, we feel that the patterns derived and dis- 
cussed in Section |] have some general value from a didactical point of view. They illustrate 
the image distortion not only in cosmological models but in all spacetimes of the respective 
Petrov type (in lowest non-trivial order with respect to the distance) and the dependence of 
this effect on the observer's velocity. This might be helpful for associating the Weyl tensor 
with some geometrical imagination. 



2 Derivation of the distortion field 

As already mentioned in the introduction, distortion as a cosmological observable was in- 
troduced in a pioneering paper by Kristian and Sachs ||. In this section it is our goal to 
rewrite the basic equations by which the distortion field is described in terms of Newman- 
Penrose coefficients. In the following section we will then discuss these results for each 
Petrov type. If the reader is not familiar with the Newman-Penrose formalism and with the 
Petrov classification he or she may consult Chandrasekhar |]23|| , Chapters 1.8 and 1.9, for a 
detailed introduction. We adopt the same sign and factor conventions as Chandrasekhar. 
In particular, we use the signature (+, — , — , — ). 

In an arbitrary spacetime (i.e., a 4-dimensional Lorentzian manifold), we fix a point 
p and a pseudo-orthonormal frame (E\, E 2 , E 3 , E4) at p. The latter may be expressed 
equivalently in terms of a null tetrad (£, n, m, rh) via 

E a x = -±(m a + m a ) , E a 2 = -^(m a - rh a ) , 



cf. Chandrasekhar [p3| |, Chapter 1.8. We interpret the timelike vector E4 as past-pointing 



so that — E4 may be viewed as the 4-velocity of an observer. Then the totality of lightlike 
geodesies issuing from p into the past is in natural one-to-one relation with the set of past- 
pointing lightlike initial vectors 

k a = sim?(cos<^ E{ + sin<^ E%) + costf E% + El 
= ^{sintffVW + e- iv rh a ) + ( 1 + cos$)£ a + ( 1 - costf)n a } , (2) 

where $ and ip have their usual range as standard coordinates on the 2-sphere. We may 
interpret $ and ip as coordinatizing the celestial sphere of the observer with 4-velocity — E4 
at p. Every point on the past light cone of p can be written as exp(sfc) with some parameter 
s > and some k given by (0); if we restrict to sufficiently small s, this representation is 
unique. Henceforth we shall use the affine parameter s as a measure for the distance from 
p, and it is our goal to discuss the distortion effect in lowest non-trivial order with respect 
to s. To within this approximation, s can be replaced with the angular diameter distance 
D = s + 0(s 2 ), with the luminosity distance D = s + 0(s 2 ), or with any other reasonable 
measure of distance used in cosmology. The redshift z is related to s by an equation of 
the form z = Hs + 0(s 2 ) with a 'Hubble constant' H which, in general, is a non-constant 
function on the sky. Also, it is important to realize that z depends not only on the velocity 
of the observer but also on the velocity of the light source. (Writing cosmological relations 
in terms of power series was a basic idea introduced by Kristian and Sachs . A discussion 
of the relations between s, D, D and z needed here can be found, e.g., in Hasse and Perlick 

m.) 
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At each point of the celestial sphere, the two orthonormal vectors 



dk a 

— = i { cos ^( e ^m a + e- l *m a ) - sim?(f - n a )\ , 

span the tangent space to the celestial sphere. (Here and in the following it goes without 
saying that one has to mind the familiar coordinate singularities at $ = and i? = it.) 
Hence, every tangent vector to the celestial sphere is of the form 

dk a 1 dk a 

t = u— + v 



d$ sini? d(f (4) 
= -L { [u costi + iv) e iip m a + (u cos$ - iv) e~ iLp rh a - u sin$ (t - n a ) } 

with u, v E R. Restricting w and f to a circle, w 2 + t> 2 < r 2 , gives an 'infinitesimally thin' 
bundle of light rays with an initially circular cross-section. When following those rays into 
the past we will observe, in general, that the cross-section does not stay circular but becomes 
elliptical, owing to the shear produced in the bundle by the spacetime geometry. (Please 
note that, according to general relativity, the question of whether an 'infinitesimally thin' 
light bundle has a circular shape has an observer-independent answer; this was demonstrated 
already in 1932 by Kermack, McCrea and Whittaker f25|, cf. Schneider, Ehlers and Falco 
p6| , Section 3.4.1.) Correspondingly, if the cross-section of the bundle is supposed to be 
circular at some affine distance s, then we have to restrict u and v to a certain ellipse. 
The equations necessary to calculate the shape of this ellipse can be found in the paper 
by Kristian and Sachs ||. Their analysis is based on the well-known fact, established in 
a fundamental paper by Sachs |1[], that the shear of an infinitesimally thin bundle of light 
rays is governed by the Weyl tensor along the central ray. More precisely, their result reads 
as follows. To determine an infinitesimally thin bundle which has a circular cross-section at 
some affine parameter distance s from the event p, one has to consider a certain quadratic 
form p a bt a t b on the two-dimensional space of tangent vectors (H). In lowest non-trivial order 
with respect to s, p a b is given as 

p ab = -iC acbd k c k d + 0(s 3 ), (5) 

where C ac bd denotes the Weyl tensor at the event p, cf. eq. (28) in Kristian and Sachs 0. If 
we neglect the 0(s 3 )-term, we can read from @ that p a b is (real-symmetric and) trace- free; 
so it has two real eigenvalues which differ by sign, say y£ and — ^-e with e > 0. To get a 
bundle with circular cross-section at s one has to restrict t a to an ellipse whose major axis 

2 2 

is the eigenspace of p a b to the eigenvalue — and whose eccentricity is equal to ^e. 

The special form of the 0(s 3 )-term in (|5]) will be of no interest to us since we will restrict 
to affine distances which are small enough to neglect this term. However, it is interesting 
to note that, if this term is written as a power series in s, the factor in front of s l+2 is a 
linear function of the % th covariant derivative of the Weyl tensor at p. Hence, if the Weyl 
tensor is covariantly constant, then the 0(s 3 )-term in is zero. Correspondingly, if the 
Weyl tensor varies but little this term may be neglected even for values of s that are not 
small. 

Now we restrict to light sources at a fixed affine distance s and we assume that, for this s, 
the 0(s 3 )-term in (|5|) can be neglected. Then the distortion effect is completely determined 
by the Weyl tensor at p. Knowing C ac M allows to calculate the eigenvectors and eigenvalues 
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<p - 

Figure 1: At each point of the celestial sphere, the distortion length element is determined 
by the non-negative number e and the angle x (mod it) 

of d^D, for all values of $ and ip (hidden in k c k d ). The result may be graphically represented 
by assigning to each point of the celestial sphere a 'length element' whose direction indicates 

2 2 

the eigenspace of p ab to the eigenvalue — y£ and whose length is a measure of ^-e (for fixed 
s). We find it convenient to choose the length elements proportional to y/e, see Figure 
0. Using the notation of this figure, the 'distortion length element field' is determined by 
giving the number e and the angle x f° r each point of the celestial sphere. Kristian and 
Sachs used the quantity e = 1 + 2e as a measure for the magnitude of the distortion 
effect. Penrose and Rindler ||, Volume II, Section 8, discussed the direction field spanned 
by our 'distortion length element field', and they called this the fingerprint direction field 
of the Weyl tensor. 

If we exclude the trivial case that the Weyl tensor vanishes at p, then there are exactly 
four (not necessarily distinct) points on the celestial sphere at which e is zero. These 
four points correspond to the four principal null directions of the Weyl tensor, and the 
question of how many principal null directions coincide leads to the Petrov classification, 
see Chandrasekhar f23fl , Chapter 1.9. In this paper it is our main goal to discuss the 
distortion field over the whole sky for the various Petrov types, thereby extending the results 
of Penrose and Rindler |J, Volume II, Section 8, who analyze the fingerprint direction field 
near a principal null direction. To that end it will be convenient to express the distortion 
field (i.e., e and x as functions of d and tp) with the help of the Newman-Penrose coefficients 

= -C abcd tm b tm d , = -C abcd tnH c m d , ^ 2 = -C abcd tm b m c n d , 

^3 = -C abcd tn b m c n d , ^ 4 = -C abcd n a m b n c m d . ' ' ' 

Inserting (Q) and (|4j) into (|5]) with the 0(s 3 )-term neglected, and using the well-known 
symmetry properties of the Weyl tensor, allows one to express the quadratic form p a bt a t b 
in terms of the five complex Newman-Penrose coefficients (||). This gives a rather long 
expression, but with some elementary algebra it can be conveniently rewritten as 

Pab t a t b = ^Re{(u + zv) 2 ee- 2l x} (7) 

where e > is the modulus and —2% (mod 2-k) is the argument of the complex number 

£e -2i X = I^ (i + CO s^) 2 e 2i¥ ' + 2^ 1 sin^(l + cos^)e^ + 
3 ^ 2 sin 2 tf + 2 ^ 3 sin d (1 - cos d) e~ iv + \ # 4 (1 - cos e -2 ^ . ^ 
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Eq. (0) can be rewritten in matrix notation as 



2 2 

From this equation we can easily deduce that the eigenvalues of p a b are ^£ and — ^e, and 
that the eigenspace to the negative eigenvalue is spanned by 

• dk a 1 dk a 

t =- sm *w + cosx ^w (10) 

This shows that £ and x? as defined by have indeed the same meaning as introduced 
before, see Figure [1|. In other words, if we know the five Newman- Penrose coefficients 
^o, . . . , ^4 (i.e., if we know the components of the Weyl tensor with respect to our null 
frame), then (f|) gives us the distortion field, i.e., e and \ as functions of $ and <p. 

The distortion field depends, of course, on the frame chosen. In terms of the pseudo- 
orthonormal frame (E\, E 2 , E 3 , E4), which is related to the null frame (£., n, m, fti) via ([!]), 
the behavior of the distortion pattern under Lorentz transformations can be summarized 
in the following way. (We are interested only in proper Lorentz transformations, i.e., those 
which do not change the temporal or spatial orientation.) Changing the spatial vectors 
Ei, E2, E3 and leaving the timelike vector E4 unaffected leads to a rigid rotation of the 
celestial sphere. Changing E4 leads to a conformal transformation of the celestial sphere, 
see Penrose and Rindler ||, Volume I, for a detailed discussion, and to a corresponding 
deformation of the distortion pattern. We want to discuss this dependence on the observer 
in some detail for the various Petrov types in the next section. As our analysis will be 
based on the representation of the distortion field in terms of Newman- Penrose coefficients, 
we shall need the behavior of those coefficients under proper Lorentz transformations. The 
relevant formulae are conveniently listed in the book by Chandrasekhar [p3fl , Chapter 1.8, 



who makes use of the fact that all proper Lorentz transformations can be composed of 
(spacetime) rotations of class I, II, and III, defined respectively by 



I : £ > £, m 1— > m + a£, ffi 1— >■ m + a£ , tihtiI am + am + \a\ 2 £ , (11) 
II : nun, mwm + k, m fh + bn , £ 1— ► £ + bm + bm + \b\ 2 n , (12) 
III : £ h-> A~ l £ , nH An, m^ e i0 m , rn h-> e~ i0 m , (13) 

with complex numbers a and b and real numbers A > and (mod 2n) . The transformation 
behavior of the Newman- Penrose coefficients is for rotations of class I: 



^o^^o, ^1 ^1 + a^o , ^2 ^ ^2 + 2a\l>i + a 2 \l>o , 
$ 3 h$ 3 + 3a^ 2 + 3a 2 ^i + a 3 ^ , ^4 >-> ^4 + 4a^ 3 + 6a 2 ^ 2 + 4a 3 ^i + a 4 ^ ; 

for rotations of class II: 

^0 ^ ^0 + 46^i + 66 2 ^ 2 + 46 3 ^ 3 + 6 4 ^ 4 , ^ 1 i-> ^1 + 36^ 2 + 36 2 * 3 + 6 3 ^ 4 , 
^2 ^ ^2 + 26^3 + 6 2 ^ 4 , \t> 3 i-> ^ 3 + W> A , * 4 ^ f 4 ; 



(14) 



(15) 



and for rotations of class III: 



^0 >-> A^e^^ , ^ A~ e \l>i , ^ 2 ^ *s 
^ 3 i-> Ae~ ie ^ 3 , ^ 4 i-> A 2 e" 2ie ^ 4 . 



(16) 
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3 Discussion of the distortion field 



Eq. (IS) gives the distortion field, i.e., e and x as functions of $ and ip, in an arbitrary frame. 
It is our goal to discuss (H) for each Petrov type. To that end we first choose a frame such 
that ^4 7^ which can always be achieved by a rotation of class I, see fll4"|). (Here and in 
the following we exclude the trivial case where the Weyl tensor is zero, i.e., of Petrov type 



O.) We then apply a rotation of class III such that ^4 = 2, see ([U|). In the resulting frame, 
(H) is equivalent to 

4 

£e ~ 2lX = H ( z n ^ 1 + cos ^ e%iPl2 - v 7 ! - costf e~ iLpl2 ) (17) 

where z%, z 2 , Z3, z 4 are the roots of the equation 

^ + 2* 3 * 3 + 3* 2 z 2 + 2M/ 1 *+!*o = 0. (18) 

(To prove this, set the left-hand side of (p~8| ) equal to [z — Z\){z — z 2 )(z — z 3 )(z — z±) and 
verify that this yields the same equations for z%, z 2 , z 3 , z± as equating the right-hand sides 
of (H) and (|i~7|) .) From this representation we immediately read that the distortion field 
has exactly four (not necessarily distinct) zeros, corresponding to the four principal null 
directions of the Weyl tensor. In particular, we read that at the 'north pole' $ = there is 
a zero of multiplicity k if and only if 



m i = for % = 0, . . . , k - 1 
^0. 



(19) 



Thus, if we know that there is a principal null direction of multiplicity k, then we know 
that, in an appropriately chosen frame, the distortion pattern is given by @ with (|19|) and 
^4 = 2. We shall call such a frame adapted to the principal null direction henceforth. 

Before discussing the global features of the distortion patterns for each Petrov type, we 
briefly show that our previous results easily allow one to analyze the local structure of the 
distortion pattern near a principal null direction of multiplicity k, thereby reproducing the 
results of Penrose and Rindler [Q, Volume II, Section 8. After inserting, for k = 4,3,2, 1, 
(|19"D and ^4 = 2 into (§) we may analyze the resulting pattern near $ = 0, see Figure 0. In 
each picture $ gives the distance from the center and ip is the azimuthal angle encircling the 
center; the length elements are proportional to ^Je. By inspection we find that for k — 4, 3, 1 
the pattern is uniquely determined, in a neighborhood of = 0, up to diffeomorphisms, i.e., 
the differential-topological structure of the distortion pattern is unique. For k = 2, the 
pattern depends on whether or not, in the representation of (|l~9|) , the argument of ^2 is zero 
(mod 2ir), i.e., ^2 is real and positive. If the argument of ^2 is zero (mod 2n), then the 
integral lines of the distortion field encircle the double zero, otherwise they are spiralling into 
the double zero, with a twist depending on the argument of ^2- Moreover, from Figure 
we read that the distortion field is generated by a continuous vector field near a principal 
null direction of multiplicity k if and only if k is even. In other words, only for k = 4 and 
k = 2 is it possible to assign an orientation to the distortion length elements in a consistent 
way. These features have already been discussed by Penrose and Rindler. In particular, the 
reader should compare Figure 8-3 in Penrose and Rindler , Volume II, with our Figure |2|. 

We shall now discuss the explicit and global structure of the distortion pattern for each 
Petrov type, rather than the local differential-topological structure near a principal null 
direction. For each Petrov type, we always start out with a frame adapted to a principal 
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Figure 2: Distortion pattern near a principal null direction of multiplicity k, cf. Penrose 
and Rindler |J, Volume II, Figure 8-3. This picture, like all the following ones, has been 
produced with MATHEMATICA. Choosing in all pictures the length elements proportional 
to y/e, rather than to e, was motivated by the fact that otherwise their length would vary 
too strongly, in particular in the case that there is a principal null direction of multiplicity 
4. 
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null direction of highest multiplicity, and we denote the Newman- Penrose coefficients in this 
particular frame by \l/o = 0, ^i, ^3> an d ^4 = 2. We then apply an arbitrary rotation 
of class I followed by an arbitrary rotation of class III, see ([[4]) and (|16D . Under such a 
transformation, which is the general form of a proper Lorentz transformation keeping the 
direction of £ fixed, the Newman-Penrose coefficients change according to 

^ = 0, = A" 1 e ie \i>i , * 2 = $2 + 2afi, 

$ 3 = Ae" i8 ($3 + 3af 2 + 3a 2 * 1 ), (20) 
^ 4 = A 2 e~ 2iB (2 + 4a^ 3 + 6a 2 $ 2 + 4a 3 $ 1 ), 

see (|i~4[) and ([16|) . Inserting the new Newman-Penrose coefficients into (§) and allowing A to 
run over R + , over [0, 27r[ , and a over C will give us the distortion pattern in an arbitrary 
Lorentz frame, with the only restriction that we keep the zero of highest multiplicity at 
the north pole. Clearly, this is just a restriction on the choice of the spatial axes which 
is a matter of convenience only. Similarly, the effect of the angle is trivial in the sense 
that it produces just a spatial rotation that can be compensated for by a transformation 
ip i — > ip + const. 

3.1 Type N 

Type N is characterized by the fact that all four principal null directions coincide. In a 
frame adapted to this fourfold principal null direction we have = ^2 = ^3 = (and, of 
course, if>o = and ^4 = 2). Inserting the transformed Newman-Penrose coefficients ( |20| ) 
into (H) gives us the distortion pattern in an arbitrary Lorentz frame (except for spatial 
rotations), 

e e~ 2ix = A 2 e ~ 2i{0+v) (1 - cos i9) 2 , (21) 

see Figure |3|. In the hatted frame (i.e., for A — 1 and = 0), ( pl|) involves no parameter 
whatsoever. Hence, the distortion pattern for type N is universal in the sense that at any 
two points (in the same or in different spacetimes) where the Weyl tensor is of type N we 
can choose observers that see exactly the same distortion patterns. Changing to some other 
observer has the only effect of introducing a scaling factor A which is constant over the sky. 
(The angle just produces a spatial rotation.) By choosing a sequence of observers whose 
4-velocities approach the fourfold principal null direction, the distortion effect can be made 
arbitrarily large, A — > 00 . By choosing a sequence of observers whose 4-velocities approach 
some other null direction, the distortion effect can be made arbitrarily small, A — > . 

3.2 Type III 

Type III is characterized by the fact that three principal null directions coincide whereas 
the fourth is different. In a frame adapted to the threefold principal null direction we have 
$1 = f 2 = 0, $3 ^ (and \l/o — 0, ^4 = 2). Inserting the transformed Newman-Penrose 
coefficients (p0|) into (§) yields 

ee- 2ix = 2Ae- i{e+ ^ ^ 3 sintf (1 - cos??) + A 2 e~ 2i ( e+ ^ (1 + 2a^ 3 ) (1 - cos??) 2 , (22) 

see Figure f|. Whatever the (non-zero) value of ^3 may be, by choosing a, A and ap- 
propriately we may give any two values we like to the coefficients Ae~ % ®^3 G C \ {0} and 
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Figure 3: Distortion pattern for Petrov type N. In this picture, and in the following ones, 
the celestial sphere is given in Mercator projection, with if as the horizontal and $ as the 
vertical variable. The picture shows the pattern given by (|2lD with Ae~ l@ = 0.05. 



A 2 e~ 2ie (l + 2a^ 3 ) e C in Thus, type III shows the same universality property as type 
N: For any two points where the Weyl tensor is of type III we can find observers that see 
exactly the same distortion patterns. However, the effect of a Lorentz transformation is now 
more complicated than for type N. Choosing a such that 1 + 2a \&3 = gives us a frame 
such that the simple zero of the distortion pattern is at the south pole (i? = 7r), i.e., at the 
antipodal point of the threefold zero. This restriction on the frame, which is equivalent to 
requiring that = -^(£ + n) lies in the plane spanned by the two different principal null 
directions, still allows to choose A (and G) arbitrarily; changing A has the same magnifying 
or demagnifying effect on the distortion pattern as for type N (and 6, as always, produces a 
trivial spatial rotation). If we now choose a different value for a, thereby boosting the frame 
such that -E4 is no longer in the plane spanned by the two different principal null directions, 
the simple zero of the distortion pattern moves away from the south pole. Hence, at a 
type-Ill spacetime point the observer can easily read from the distortion pattern whether 
his or her 4-velocity is in the plane spanned by the two different principal null directions. 

3.3 Type D and type II 

In a frame adapted to a principal null direction of multiplicity two the Newman-Penrose 
coefficients satisfy $1 = 0, $2 / (and \l/o = 0, ^4 = 2). The remaining two principal 
null directions are given by the two non-zero solutions of equation ( |T8| ) with the hatted 
coefficients, i.e., 

Z3/4 = -*3 ± \/^-3^2 • (23) 

If these two solutions coincide, the Weyl tensor is of type D, otherwise it is of type II. Hence, 
we have 

Type D : 3 ^ 2 = ^ (24) 
Type II : 3 ^ 2 ^ ^l- (25) 
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Figure 4: Distortion pattern for Petrov type III, seen by an observer whose 4-velocity lies in 
the plane spanned by the two different principal null directions (top) and by another observer 
(bottom). The first picture shows the pattern ee~ 2lx = 0.007 siwd ( 1 — cos$), the second 
one results by applying a rotation of class I with a = 0.52 — i 0.52, followed by a rotation 
of class III with Ae~ t@ = 0.56. Here and in the following pictures the apparently strange 
choice of numerical values is motivated by our desire to have the zeros of the distortion 
pattern on (or close to) grid points. 
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Inserting the transformed Newman- Penrose coefficients (|20| ) into (§) yields the distortion 
pattern for type D and type II in an arbitrary frame, 

£e -2i X = 3^ 2S in 2 ^ + 2Ae- 4(e+¥,) (^3 + 3a^ 2 )sin^(l-cos^) + 

A 2 e - 2l(e+Lp) (1 + 2 a ^ 3 + 3 a 2 ^ 2 ) (1 - cos d) 2 . ' ^ ' 

We first discuss type D. With fl24|) , (|26|) simplifies to 

ee" 2!X = (^ 3 sin^ + Ae- i(e+ ^(l + a^ 3 )(l-cos^)) 2 , (27) 

see Figure |^ and Figure || From this equation we read that for type D the universality 
property known from types N and III is not satisfied: The distortion pattern depends on 
the (non-zero) value of \l/3, i.e., there is a one-complex-parameter family of genuinely distinct 
type-D spacetime points. Even the differential-topological structure of the distortion pattern 
depends on ^3. If is real (i.e., if \l/ 2 = is real and positive), then the integral lines 
of the distortion field are closed curves encircling each of the two double zeros; if ^3 is 
non-real, those integral lines start at one double zero and terminate at the other. The effect 
of a Lorentz transformation is as follows. Choosing a such that 1 + a ^3 = gives us a 
frame such that the two double zeros are in antipodal positions at the sky. This is the 
case if and only if the vector E4 = -^=(£ + n) lies in the plane spanned by the two double 
principal null directions. Restricting the frames in this way leaves the freedom of choosing A 
and G arbitrarily; this, however, has no effect on the distortion pattern. In other words, all 
observers whose 4- velocity is in the plane spanned by the two double principal null directions 
see exactly the same pattern. Changing a, i.e., boosting to a frame such that E4 is not in 
the plane spanned by the two double principal null directions, has the effect that the second 
double zero moves away from the south pole. Hence, by observing the distortion pattern 
an observer can find out whether the Weyl tensor at the point of observation is of type D, 
and whether his or her 4-velocity is in the plane spanned by the two double principal null 
directions. Moreover, the differential-topological structure of the pattern alone allows to 
find out whether in a frame adapted to one of the two double principal null directions the 
Newman- Penrose coefficient \l/ 2 — §^3 is real and positive. 

We now turn to type II, see Figure |7| and Figure || By choosing a such that ^ 3 + 3 a ^ 2 = 
and choosing A and such that A 2 e~ 2lB (3 \l/ 2 — \1/|) = 9\l/ 2 we may boost to a frame 
such that E4 = -^(£ + n) lies in the plane TL spanned by the two simple principal null 
directions and the double principal null direction is spanned by the sum of E4 and a vector 
perpendicular to 7i. These properties fix E4 uniquely. In this frame, the distortion field is 
given by 

ee - 2i * = 3^ 2 (sin 2 tf + e" 2 ^ (1 - cos??) 2 ) . (28) 

In this representation the two simple principal null directions are situated at the equator 
•& = 7r/2 and opposite to each other. Similarly to the type-D case, there is a one-complex- 
parameter family of genuinely distinct type-II spacetime points, depending on the value 
of \l/ 2 . Again, the differential-topological structure of the distortion pattern depends on 
whether or not ^ 2 is real and positive. If this is true, then all the integral lines of the distor- 
tion field are closed curves encircling the double zero, with the exception of one particular 
integral line that connects the two simple zeros. If \l/ 2 is not real and positive, then every 
integral line connects the double zero with one of the two simple zeros. 
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Figure 5: Distortion pattern for Petrov type D with ^ 2 = real and positive, seen by 
an observer whose 4-velocity lies in the plane spanned by the two different principal null 
directions (top) and by another observer (bottom). We have chosen = 0.08. The second 
picture results from the first one by applying a rotation of class I followed by a rotation of 
class III with aAe~ ie = 0.58. 
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Figure 6: Distortion pattern for Petrov type D with ^ 2 — §^3 not real and positive, seen 
by an observer whose 4-velocity lies in the plane spanned by the two different principal null 
directions (top) and by another observer (bottom). This time we have chosen ^ 3 = 0.08e l7r//4 . 
As in Figure [5|, the second picture results from the first one by applying a rotation of class 
I followed by a rotation of class III with aAe"' e = 0.58. 
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Figure 7: Distortion pattern for Petrov type II with \& 2 real and positive, seen by an 
observer with the special 4-velocity such that the distortion pattern is given by ( ^8l) (top) 
and by another observer (bottom). We have chosen \& 2 — 0.001. ^3 is arbitrary except for 
3^|/ 2 -L The second picture results from the first one by applying a rotation of class I 
with a = 0.32 + i 0.73, followed by a rotation of class III with Ae~ l@ = 0.93 + i 0.33. 
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Figure 8: Distortion pattern for Petrov type II with ^ 2 not real and positive, seen by an 
observer with the special 4-velocity such that the distortion pattern is given by ( p8"D (top) 
and by another observer (bottom). This time we have chosen ^ 2 = O.OOle"' 7 ' 4 . Again, 
^3 is arbitrary, except for 7^ and the second picture follows from the first one by 
applying a rotation of class I with a = 0.32 + i 0.73, followed by a rotation of class III with 
Ae~ ie = 0.93 + i 0.33. 
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3.4 Type I 



Type I is characterized by the fact that all four principal null directions are distinct. In a 
frame adapted to one of them we have ^ (and \l/o = 0, ^4 = 2). With the transformed 
Newman-Penrose coefficients ( p0|) we find the following expression for the distortion pattern 

ee - 2ix = 2A-V (0+v) ^isintf (1 + costf) +3(^ 2 + 2a^i)sin 2 ^ + 

2 A e' l{e+tp) (# 3 + 3 a ^ 2 + 3 a 2 ifj) sin d (1 - cos •&) + (29) 
A 2 e ~ 2i{@+Lp) (1 + 2 a ^3 + 3 a 2 ^ 2 + 2 a 3 #1) (1 - cos ) 2 . 

We may choose a such that 1 + 2 a ^3 + 3 a 2 ^2 + 2 a 3 ^1 = 0, thereby moving the second 
zero to the south pole 1? = it. With a fixed that way, we may choose A and G such that 
A~ 2 e 2l@ ^i = (^3 + 3 a ^2 + 3a 2 \I/i), thereby moving the remaining two simple zeros in 
symmetric positions, $3 + $4 = 7r and ^3 + y?4 — 27r. The resulting frame is characterized 
by the fact that E± = -^=(£ + n) lies in the plane TL spanned by two of the simple principal 
null directions and the other two simple principal null directions are spanned by vectors of 
the form E4 + X + Y and E4 + X — Y where X and Y are perpendicular to E4 and X 
is perpendicular to 7i. These properties characterize E4 uniquely, except for the possibility 
of permutating the four principal null directions. In this frame, the distortion pattern takes 
the form 

£ e~ 2ix — Q sin 1? ((1 + cos 1?) e itp + (1 — cost?) e~ iip ) + Psin¥ (30) 

where Q = 2A~ 1 e i@1 & 1 G C \ {0} and P = 3 ( ^ 2 + 2a$ 1 ) G C are parameters that 
characterize the spacetime point. Thus, there is a two-complex-parameter family of type-I 
spacetime points which are genuinely distinct in the sense that their distortion patterns are 
not related by a Lorentz transformation (i.e., by a conformal transformation of the celestial 
sphere). However, the distortion patterns of type-I spacetime points are always related by 
diffeomorphism, i.e., the differential-topological properties of the distortion pattern are the 
same for all values of Q and P. The top part of Figure ^| shows the distortion pattern given 
by fl30l ) for some values of Q and P; the bottom part demonstrates the effect of a Lorentz 
transformation on this pattern. Two real parameters are necessary for fixing the positions 
of the zeros in the first picture. The other two real parameters in Q and P fix the overall 
scaling factor and the twist of the pattern. 

As an aside, we mention that there is always a frame in which the four simple zeros on 
the sky form a disphenoid, i.e., a tetrahedron such that opposite edges have equal length, 
see Penrose and Rindler ||, Volume II, Section 8. However, with our convention of keeping 
one of the zeros at the north pole the distortion field in this particular frame is given by a 
rather awkward expression. 

4 Conclusions 

In this paper we have given a new formula, eq. @, for the distortion effect on light sources 
around an observer in an arbitrary spacetime. This formula, which is exact to within lowest 
non-trivial order with respect to the distance between light sources and observer, comes 
about by rewriting classical results of Kristian and Sachs in terms of the Newman- 
Penrose formalism. On the basis of this formula we have discussed the distortion patterns 
graphically for each Petrov type. We feel that the resulting pictures have some didactic value 
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Figure 9: Distortion pattern for Petrov type I, seen by an observer with the special 4- 
velocity such that the distortion pattern is given by (|30"D (top) and by another observer 
(bottom). We have chosen P = -0.0018 -i 0.0056 and Q = -0.0022 -i 0.0035. The second 
picture results from the first one by applying a rotation of class I with a = —0.25 + i 0.14, 
followed by a rotation of class III with Ae~ l@ = —1.37 + i 0.41. 
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as visualizing the effect of the Weyl tensor on light rays. In this respect our results extend 
the work by Penrose and Rindler H] on the 'fingerprint direction field' of the Weyl tensor 
insofar as (i) we discuss not only the direction but also the magnitude of the distortion, 
(ii) our patterns are metrically correct on the whole celestial sphere, not just differential- 
topologically near a principal null direction, (iii) we discuss the dependence on the observer's 
velocity. 

As a special application, it is an interesting question to ask whether the distortion 
patterns derived in this paper have some relevance in view of cosmology. As mentioned 
already in the introduction, image distortion produced by large-scale structure has been 
observed recently [pi], [Tl], 0, |13[ [yj by statistically evaluating galaxy ellipticities in selected 
fields on the sky. Following the original ideas of Kristian and Sachs 0, we want to ask 
whether observations of this kind can be used to determine the Weyl tensor 'here and 
now' in an appropriately smoothed universe. For discussing this question we subdivide our 
celestial sphere into fields of some specified size and we specify a maximal distance up to 
which galaxies are to be observed. By averaging galaxy ellipticities over each field and 
assigning the result to the field center we get the distortion pattern of a smoothed universe. 
Here the smoothing procedure refers to averaging the Weyl tensor over conic regions, with 
the observer at the tip of the cone, which is not directly related to averaging the matter 
density. If the field size and/or the maximal distance has been chosen small, then the 
distortion effect in each conic region is essentially determined by local inhomogeneities. 
In other words, even in the smoothed universe the Weyl tensor will have a fairly strong 
variation, so the 0(s 3 )-terms in (|5]), which involve covariant derivatives of the Weyl tensor, 
cannot be neglected. For this reason we cannot expect that the resulting pattern coincides 
with one of the patterns derived in the preceding section. However, it does make sense to 
compare the observed distortion pattern with the ones derived in the preceding section if we 
increase the field size and the maximal distance sufficiently such that the distortion effect 
of local inhomogeneities becomes irrelevant. There are three possible outcomes. 

First, the distortion effect may vanish over the whole sky, for field size and maximal 
distance chosen appropriately. This would indicate that, at this level of averaging, the Weyl 
tensor is zero. This is what we expect if we rely on the assumption that our universe differs 
from a Robert son- Walker cosmos only by local 'lumps' whose effect on light rays averages 
to zero. The present observations are in accordance with this assumption, but one has to 
keep in mind that they cover only a very small portion of the sky. When supported by 
future observations that cover large parts of the sky, this can be viewed as an independent 
confirmation of the standard cosmology assumptions. It is remarkable that this line of 
thought is based on geometry only, i.e., Einstein's field equation is not used. 

Second, the distortion pattern may become similar to one of the patterns derived in the 
preceding section, for some field size and some maximal distance. This would indicate that 
in our standard cosmology model the Robertson- Walker background has to be replaced by 
some other background model, namely by one with a non-trivial Petrov type. In this case 
we could read from the distortion pattern not only the Petrov type of the Weyl tensor 'here 
and now' (in the smoothed universe) but also the position of our 4- velocity vector in relation 
to the principal null directions. 

Third, the distortion effect may be non-zero but different from all the patterns derived in 
the preceding section, for any field size and any maximal distance. In this case there would 
be the following possible explanations, (i) The galaxies inside one field have non-randomly 
distributed actual axes, even for large field size and large maximal distance. This could 
be taken as indicating a prefered direction in the universe such as given, e.g., by a cosmic 
rotation or by a cosmic magnetic field, (ii) Even in the smoothed universe the derivative 
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terms of the Weyl tensor are so large that the 0(s 3 ) terms in @ cannot be neglected, (iii) 
In some fields the average distance of the galaxies is larger than in other fields. As long 
as the 0(s 3 ) terms can be neglected, this has an effect only on the magnitude but not on 
the direction of the distortion effect, i.e., it has an effect on e, as a function on the celestial 
sphere, but not on \. 

In this sense, observing the distortion effect can provide us with some information on 
the Weyl tensor in the universe, based on geometry arguments alone. This line of thought is 
meant as complementary to the ongoing efforts of using distortion observations for tracing 
the matter distribution in the universe, based on the approximative assumptions of the 
weak-lensing formalism. 

Two questions are left open for future work. The first one is whether there are some 
classes of spacetimes, e.g. Bianchi models, for which the distortion patterns can be cal- 
culated exactly, not just to within lowest order with respect to the distance. The second 
one is whether distortion patterns can be linked to the matter distribution without making 
approximative assumptions of the weak-lensing kind. It is true that, quite generally, the 
Weyl tensor is that part of the curvature which is not determined, at a particular event, 
by the energy-momentum tensor at that event via Einstein's field equation. However, the 
divergence of the Weyl tensor is related to derivatives of the energy-momentum tensor by 
the equation C abcd -d = T c t a;b l — lg c [ a T ; ^ which follows from the Bianchi identity and Ein- 
stein's field equation (see Kundt and Triimper |27j or Ellis 01). Maybe this equation can 



be used, at least in special classes of spacetimes, for gaining some information about the 
energy-momentum tensor from distortion patterns. In this connection recent work by Frit- 
telli, Kling and Newman [^9], is of some interest. In these articles the authors discuss 
image distortion in arbitrary spacetimes without using approximative assumptions of the 
weak-lensing kind and also without using series expansions of the Kristian-Sachs type. 
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